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We present a large family of exact solitary wave solutions of the one dimensional Gross- 
Pitaevskii equation, with time- varying scattering length and gain/loss, in both expulsive and regular 
parabolic confinement regimes. The consistency condition governing the soliton profiles is shown to 
map on to a linear Schrodinger eigenvalue problem, thereby enabling one to find analytically the 
effect of a wide variety of temporal variations in the control parameters, which are experimentally 
realizable. Corresponding to each solvable quantum mechanical system, one can identify a soliton 
configuration. These include soliton trains in close analogy to experimental observations of Strecker 
et al, [Nature 417, 150 (2002)], spatio-temporal dynamics, solitons undergoing rapid amplification, 
collapse and revival of condensates and analytical expression of two- soliton bound states, to name 
a few. 
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Coherent atom optics is the subject of much current 
interest due to its relevance to both fundamental aspects 
of physics, as well as to technology 0, 0, 0, 0, 0, EJ • For 
that purpose, lower dimensional condensates e.g., cigar- 
shaped Bose-Einstein condensates (BECs) have been the 
subject of active study in the last few years m. l& M. [ lCI 
Hlj . Observations of dark and bright solitons |l2lll3l [l3|. 
particularly the latter one, since the same is a conden- 
sate itself, have generated considerable interest in this 
area. This has spurred intense investigations about the 
behavior of condensates in the presence of time-varying 
control parameters. These include nonlinear ity, achiev- 
able through Feshbach resonance 

gain/loss 

and the oscillator frequency The fact that for a 

condensate in oscillator potential, exact solutions of the 
Gross-Pitaevskii (GP) equation are not available, makes 
it extremely difficult to examine the effects of time varia- 
tion in the aforementioned parameters. In the context of 
pulse propagation in non-linear optical fibers, a number 
of authors have recently investigated the effects of vari- 
able non-linearity, dispersion and gain or loss. Moores 
analyzed this problem for constant dispersion and non- 
linearity and a distributed gain |l9j |. whereas Kruglov et 
al, considered the same problem with all the param- 
eters in a variable form 20]. Serkin et al, have de- 
rived the nonlinear Schrodinger equation (NLSE) with 
distributed coefficients as compatibility condition of two 
first order equations, demonstrating the applicability of 
inverse scattering transform method to this type of prob- 



* atreOprl . res . in 
tprasanta@prl . res . in 
■i-girish. agarwal@okstate . edu 



lems 21]. They write down the general equation relating 
the distributed coefficients with the solution parameters 
and obtain a number of exact solutions. Recently, exact 
solutions of a driven NLSE with distributed dispersion, 
nonlinearity and gain, which exhibits pulse compression 
in a twin-core optical fiber, has also been obtained [22^ . 

In this Letter, we present a large family of exact solu- 
tions of the quasi one-dimensional GP equation, which is 
the familiar NLSE, with time varying scattering length, 
gain/loss, in the presence of an oscillator potential, which 
can be both expulsive or regular. It is shown that, 
the consistency condition governing the soliton profiles 
identically maps on to the linear Schrodinger eigenvalue 
problem, thereby allowing one to solve analytically the 
GP equation for a wide variety of temporal variations in 
the control parameters. Corresponding to each solvable 
quantum mechanical eigenvalue problem, one can iden- 
tify a soliton-like profile. These can be dark or bright 
and the oscillator frequencies can have a variety of tem- 
poral profiles, including linear, quadratic, exponential, 
smooth step-like potential and kicked oscillator scenario. 
Our solutions exhibit soliton trains, spatio-temporal dy- 
namics of solitons, the formation of two-soliton bound 
states, to name a few. Interestingly, it was found that 
some of these solitons can periodically exchange atoms 
with a background. Further, our analytical results are 
closely related to experiments [3 Q . Amplification of 
atomic condensate and condensate compression are ob- 
served in the parameter domain, which are amenable for 
experimental verification. 

Analogous to the experimental observations of Strecker 
et al, we obtain bright soliton trains in the presence of 
harmonic confinement [3] • Some of the soliton solutions 
in the presence of regular harmonic confinement, with 
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appropriately tailored gain/loss, exhibit collapse and re- 
vival phenomena, with an increase in amplitude. We 
treat the attractive regime rather exhaustively since the 
bright solitons are themselves condensates. The analytic 
expression for the dark soliton in the repulsive sector has 
been presented, detailed analysis of the same can be car- 
ried out in a straightforward manner. 

We start with a zero temperature BEC of atoms, con- 
fined in cylindrical harmonic trap Vo(x,y) = muo\(x 2 + 
y 2 )/2 and a time-dependent harmonic confinement, 
which can be both attractive and expulsive, along the 
z-direction V\(z,t) = m(J^(i)z 2 /2: 



ih 



d V(r,t) 
dt 



2m 



V 2 + Z7|^(r,£)| 2 + F + i 



*(r,t), 
(1) 

where U = Anh 2 a s (t) /m and V = Vb+Vi. It is worth not- 
ing that, condensate can interact with the normal atomic 
cloud through three-body interaction, which can be phe- 
nomenologically incorporated by a gain/loss term rj(t). 

To reduce Eq. JIJ to an effective one dimensional equa- 
tion, we assume that, the interaction energy of atoms is 
much less than the kinetic energy in the transverse direc- 
tion El: 



*(r,*) 



x exp —icj±t 



y 



2a\ 



(2) 



In dimensionless units, GP equation then reduces to the 
following one dimensional nonlinear Schrodinger equa- 
tion: 



1 



d zz ^j + 7 (tM 2 ^- 



1m(^V + <^-(3) 



Here, 7 = 2a a (t)/a B , M(t) = u 2 (t)/ul, g(t) = 
rj(t)/fuj±, a± = (H/muj^) 1 / 2 and a B is the Bohr's ra- 
dius. For the sake of generality we have kept M(t) time 
dependent, a constant M(t) implies an oscillator poten- 
tial which can be confining or expulsive for M > or 
M < 0, respectively. 

In order to discuss the cases of a variety of experi- 
mentally achievable profiles of the time dependent trap- 
ping potential and to obtain corresponding analytical so- 
lutions, we assume the following ansatz solution: 



^{z, t) = </Atf)F[A(t) {z - l(t)}\ exp[i$(z, t) + G(t)/2], 

(4) 

where, G(t) = f*g(t')dt', l(t) = f*v(t')dt' and we as- 
sume that phase has a quadratic form: 



=a{t) +b(t)z- -c(t)z 2 . 



(5) 



In the above equation, a(t) = cio + J* A 2 (t')dt', and 
c(t) is determined by a Riccati type equation: 



c t -c 2 (t) =M(t). 



(6) 



Interestingly, this equation can be expressed as a 
Schrodinger eigenvalue problem via a change of variable, 

c(t) = - rfh y>] : 



ip"(t)-M(t)ip(t)=0. 



(7) 



Taking advantage of this connection, below we show that, 
corresponding to each solvable quantum mechanical sys- 
tem, one can identify a soliton configuration. The fact 
that, Schrodinger equation can be exactly solved for a 
variety of V(i), gives us freedom to control the dynam- 
ics of BEC in a number of analytically tractable ways. 
This is one of the main results of this Letter. Although 
a host of time dependent oscillator frequencies can be 
addressed, we will only demonstrate in the text, a few 
experimentally realizable examples. 

We also find the following consistency conditions 



7(*) = 7oe- G(t)/2 ^)Mo, b(t) = A{t) \ 
A(t) = A exp 

*M-c(t)l(t)=b(t) 



(8) 



Now substitution of the ansatz Eq. (@J) in Eq. J3J) and 
using the consistency conditions Eq.®, we obtain the 
differential equation for the function F in terms of the 
new variable T = A(t)[z - l(t)]: 



F"(T) - \F(T) + 2kF 3 (T) = 0, where k = 



7o_ 
~A ' 



(9) 



We note that, the Eq. (0) possesses a variety of solu- 
tions in the form of twelve Jacobian elliptic functions e.g., 
cn(T, m), sn(T, m) and dn(T, m) etc. Here m is the mod- 
ulus parameter taking value in the interval < m < 1. 
These functions interpolate between the trigonometric 
and hyperbolic functions in the limiting cases m = and 
1, respectively. In the appropriate parameter regimes, 
a wide variety of solutions emerge which include soliton 
trains, akin to experimental observations of Strecker et 
al, and localized bright and dark solitons. 

The bright soliton trains of the form i/j(z,t) = 
y/A(t) cn [T/r , m] exp[z$0, t) + G(t)/2] exist for n > 0, 
where Tq = — Ao(m-\-m 2 ) /2jo and A = (m 2 + m — 1)/2tq. 
Here, all the coefficients are determined from Eqs. (O 
IHJ). It should be pointed out that, in the repulsive do- 
main 70 > the localized dark soliton is of the form 

~ A(t)z~ 



i/>(z,t) = y/A(t) tanh 



exp[i$(z,t)+G(t)/2], since 



the same vanishes at the origin. The cnoidal wave solu- 
tion for this is of sn(T, m) type. 

Below, we examine the formation and dynamics of non- 
linear excitations in the presence of an oscillator poten- 
tial, with a variety of experimentally achievable temporal 
modulations in the frequency. 




FIG. 1: Plots depicting the snapshots of soliton trains at times 
~ 5ms, ~ 70ms and ~ 150ms with parameters M — .25, 
7 — —0.5, to — .9999, Ao — 0.5 and modulus parameter 
m = .99999. 




FIG. 2: Figure depicting collapse and revival of atomic con- 
densate and amplification through periodic exchange of atoms 
with the background, for the parameter values M — 0.0025, 
to = 2.5, Ao = 3.125, 70 = —0.5 as in the Ref. [lj and gain 
parameters a\ — .01, a 2 — 8.25 and S = 1.5. 



Soliton trains in a confining oscillator. — Inspired by the 
experiments of Strecker et. al, [l^ in the presence of con- 
fining oscillator, with no gain/loss we obtain following 
moving soliton trains with velocity v(t) = Ao cos(Mq£): 



aM) sec(Mot) 
~A sec(M t)(z 



x cn 



^sin[M *]) 



rn 



s{l0) 



As shown in Fig. JJJ initially at « 5ms there is only 
one profile and as time progresses it breaks up in to many 
profiles giving rise to a soliton train. 
Collapse and revival of the condensate. — Inspired by 
the experiments of Strecker et al, we tailor the gain 
profile in the same range of the parameters, where 
one observes dramatic collapse and revival of conden- 
sates with an increase in the amplitude. Specifically we 
take, M(t) = Mq, g(t) = ait — a 2 s'm(dt) and <y(t) = 
aif2 ?f cos(dt) sec(M £), for which we obtain 
Aq sec(Mot). The complete solution 



7o exp L ^— -j 

from Eq. ©, A(t) 
can be written as 



(70 < 0) the effect of a step like change in the oscil- 
lator frequency which is well mimicked by the function 
M(t) = -4 + §[1 + tanh(£/2)], for which the coupling 



has the form j(t) 



7o 



(l+e 



70 < 0. As time pro- 



gresses, the coupling strength changes exponentially and 
the condensate wave function starts building up rapidly. 
In this expulsive case, like the experimentally observed 
bright soliton (T^, one obtains, 



A e 2t z 
(1 + e*)r 



exp[i$(z,t)]. (12) 



New solutions in the case of constant coupling. — A num- 
ber of interesting condensate profiles emerge in the con- 
stant attractive coupling regime, depending on the na- 
ture of the external potential. For the regular oscillator 
confinement one obtains a spectacular spatio-temporal 
pattern given by Eq. (|T3|) in the amplitude of the order 
parameter, as seen in Fig. The extreme increase in 
amplitude happens due to the presence of the sec (Mot) 
in the amplitude. 



ip(z,t) = v Aq sec(Mot)sech 



Ao sec(M t) 



to 



exp[G/2]e^ 



ip(z,t) = yj Ao sec(Mo£)sech 



Ao sec(M t)z 



TO 



exp [G(t)/2] e { 
(11) 

This is depicted in Fig.Q. In this case, sinusoidal nature 
of gain function implies a periodic exchange of atoms with 
the background. In the attractive domain, presence of a 
background surrounding the condensate has been seen in 
the experiments leading to formation of bright solitons 
[~L3L H^ . In this light, the collapse and revival of the 
condensate, having a sinusoidal exchange of atoms with 
the background, may be amenable for verification. 
Soliton formation from step changes in the trap 
potential. — Below, we explicate in the attractive domain 



(13) 

. New bound state for solitons. — Interestingly, as shown 
in Fig. (|2J), in the expulsive domain one obtains bound 
states of solitons in the same parameter regime. The an- 
alytical expression for this remains same as in Ea. (|T3| . 
except for the value of A(t) which is now given by A(t) = 
Aosech[Moi\. It is worth mentioning that, in optical 
fibers with a variable dispersion similar structures have 
been seen through numerical investigations [2^ . 

Keeping in mind, the nontrivial nature of the phase 
and its connection with soliton profile, the nature of the 
spatio-temporal dynamics of the same for various config- 
urations is worth discussing. The quadratic nature of the 
phase with respect to the space coordinate introduces a 
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FIG. 3: Figure depicting spatio-temporal pattern of conden- 
sate wave function in a confining oscillator potential (M(£) = 
0.25) with constant attractive coupling j(t) = —0.5 and 
A 0.5, to 1. 



U(.z,t) 




FIG. 4: Formation of two-soliton bound state in an expulsive 
potential( M(t) — —0.25) with all the other parameters having 
same values as in Fig.©. 

spatial chirping effect. Furthermore, its coefficient c(t) is 
also time-dependent and is connected with the external 
potential through Eq. (0; hence the trap potential has 
a signature on the same. The purely temporal part a(t) 
is also sensitive to trap potential as depicted by Eq. JSJ. 

In conclusion, for one dimensional GP equation in an 
oscillator potential, with time dependent coupling and 
gain/loss, we have obtained a wide class of exact solu- 
tions. These include experimentally observed bright soli- 
ton profile in the attractive coupling regime for the expul- 
sive case. The fact that the equation governing the soli- 
ton profiles in the presence of time dependent harmonic 
oscillator, maps to the Schrodinger equation, opens up 
a host of opportunities. Corresponding to any solvable 
quantum mechanical potential, a time dependent cigar- 
shaped BEC profile can be obtained. We have explicated 
amplification of BEC profile through smooth variation 
of the oscillator frequency. Formation of soliton bound 
states and spectacular spatio-temporal patterns that can 



manifest in this nonlinear system with time-dependent 
control parameters are demonstrated. One observes dra- 
matic compression and localization of broad condensate 
profiles, which may have technological implications. We 
have analyzed the nonlinear excitations in the presence of 
oscillator potential having linear, quadratic, Morse and 
other type of time dependencies. 
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